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This chapter reviews several quantities that have been proposed in scattering theory to charac- 
terize the temporal aspects of one dimensional collisions: the dwell time, the delay time, the decay 
time, and times characterizing transient effects or the attainment of stationary conditions. Some 
aspects of tunnelling times are also discussed. 



I. INTRODUCTION 

T-H ! 

Quantum scattering theory deals with collisions, namely, interactions which are essentially localized in time and 

■ space. This means that the interaction potential must vanish rapidly enough in coordinate space, so that the wave 
packet tends to free-motion incoming and outgoing asymptotic states before and after the interaction is effective. The 
scope of scattering theory also includes "half-collisions" or "decay processes" where the stage before the collision is 
ignored, i.e., the evolution of the system is only considered from the interaction region. 

This chapter reviews various quantities that have been proposed in scattering theory to characterize the temporal 
aspects of the collision. A quantum wave packet collision with a potential barrier in one dimension (ID) is fully 
described by the evolution of the wave function ip(x, t) from the incoming to the outgoing asymptotic states. However, 
the whole information contained in ijj(x, t) is hardly required. A few well chosen quantities are often enough to provide 
a fair picture of the dynamics. In particular, one of these elementary parameters is the transmission probability Pt, 
but to describe the time dependence we also need to quantify the duration of the collision, the arrival time at a 
' detector, the decay time of an unstable state, the asymptotic behaviour at short and large times, or response times, 
such as the time required to "charge" a well or achieve stationary conditions when a source is turned on. 

In spite of the inherent time dependence of collisions, the treatises on quantum mechanics or scattering theory 
concentrate on solutions of the time-independent Schrodinger equation. This is in part because many scattering 
experiments to obtain cross sections are performed in quasi-stationary conditions, and also because the stationary 
scattering states form a basis to analyze the actual time dependent collision. In many cases wave packet scattering 
is relegated to justify the cavalier obtention by stationary methods of cross section expressions, and occasionally to 
i-Q ' discuss resonance lifetimes. Another widespread limitation is the exclusive interest in the final results of the collision 
at asymptotic distances and times, which has been generally justified because "the midst of the collision cannot be 
observed" . However, while it is true that in many collision experiments only the asymptotic results are observed, 
modern experiments with femptosecond laser pulses or other techniques known as "spectroscopy of the transition 
state" do probe the structure and the evolution of the collision complex Q. Also, in quantum kinetic theory of gases, 
q-i accurate treatments must abandon the "completed collision" approximation and use a non-asymptotic description, e.g. 

■ in terms of Moller wave operators instead of S matrices, as in the Waldmann-Snider equation and its generalizations 
for moderately dense gases 0. 

The theory has to adapt to these new trends by paying more attention to the temporal description of the collisions. 
Even if we restrict ourselves to asymptotic aspects, the cross section does not contain the whole information available 
in a scattering process, since it is only proportional to the modulus of the ^-matrix elements. Information on the 
phase is available from delay times with respect to free motion. In fact, the full collision and not just the asymptotic 
regimes should be understood to control or modify the products. This has motivated a recent trend of theoretical 
and experimental work to investigate the details of the interaction region and transient phenomena. 

In this chapter we restrict ourselves to one dimensional scattering. Many physical systems can be described in one 
dimension: the application of the effective mass approximation to layered semiconductor structures leads to effective 
one dimensional systems ||; some surface phenomena are described by ID models and chemical reactions can 
in certain conditions be modelled by effective one dimensional potentials g. The simplicity of ID models have 
made them valuable as pedagogical and research tools. They facilitate testing hypothesis, new ideas, approximation 
methods and theories without unnecessary and costly complications. For the same reasons they are frequently used to 
examine fundamental questions of quantum mechanics. In particular, the time quantities treated in this book, such 
as tunneling or arrival times, have in most cases been examined in one dimensional models. Many results for ID are 
inspired by results previously obtained in 3D, although the direct translation is not always trivial or possible. 3D 
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collisions with spherically symmetric potentials are described, by decomposition into partial waves, on the half line, 
whereas ID collisions involve the full line and a doubly degenerate spectrum. 

The chapter is organized as follows: Section 2 provides a minimal overview of formal ID scattering theory. The 
treatment is "formal" because no mathematically rigorous proofs are given. Instead, we summarize the operator 
structure of the theory and the results needed to define characteristic times later on. For a more rigorous mathematical 
presentation see e.g. ||. Sections 3, 4 and 5 are devoted, respectively, to the dwell time, the delay time, and decay 
times (the exponential decay and its deviations). Quantities related to the tunneling time conundrum are scattered 
in several parts of the book. In particular, sections 4.1 and 4.6 discuss the Hartman effect and negative delays, while 
section 6 discusses the role of the Buttiker-Landauer "traversal time" in the time dependence of evanescent waves. A 
detailed discussion of the arrival times is left for Chapter 



II. SCATTERING THEORY IN ID 



A. Basic premises and notation 



Let H = Hq + V be the Hamiltonian operator for a single particle in one dimension, where 

m - 1 (i) 

is the kinetic energy operator in terms of the momentum operator p op ^\ and V is a "local" potential operator with 
coordinate representation 

(x\V\x') = 6[x-x')V(x) . (2) 

V(x) must vanish for large values of \x\ so that the Moller operators, defined below, exist. This may certainly be 
accomplished by finite range potentials, but spatial decays with infinite tails are also possible. 
The plane waves \p) , with coordinate representation given by 

(x\ P ) = h- 1/2 e ixp / % , (3) 

are improper eigenstatesQ of p op and Hq, normalized according to Dirac's delta function, 

(p\p')=5(p-p'). (4) 

Closure relations (or resolutions of the unit operator l op ) may therefore be written in momentum or coordinate 
representation as 

/>oo />oo 

l op = / dx\x){x\= I dp\p)(p\. (5) 



B. Basic abstract and parameterized operators. 

The state vector of the particle at time t is denoted as \ip(t)) or simply as ip(t). We shall only deal with potentials 
such that at large times in the past and future certain states ip, the scattering states, tend (in a strong sense) to freely 
moving asymtotic states ^>i n and </> ou t respectively, 

^(*) <f>in(t), t->-00, (6) 

^out(i), t^oo. (7) 



1 The subscript "op", meaning "operator", is not used for all quantum operators, but only when confussion is possible with 
ordinary numbers or functions. 
2 i.e., not in the Hilbert space of square integrable states. 
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The central objects in scattering theory are the abstract Moller operators. They link the asymptotic states with ip, 



i>(t) = n + <f> in (t) , (8) 

i>(t) = n_<t> out (t) . (9) 

Another important operator is 

s = nln+ , (io) 

which links the two asymptotes, 

0out(f) = Sftn(*)> ( H ) 

It is also convenient to introduce the auxiliary "transition" operators T± as 

r± = vn± . (12) 

The explicit definition of the Moller operators is given by infinite-time (strong) limits, 

fi± = lim e iHt/h e -iH t/h r 13 s 

t — >^oo 

The domain of these operators is the Hilbert space of square integrable states, although it is very useful to consider 
an extension that can be applied on plane waves and allows to work in a momentum representation. To this end let 
us first define the parameterized operators 

n(z) = l + G (z)T(z) , (14) 
T{z) = V + VG{z)V , (15) 

where z is a complex variable with dimensions of energy, and G(z) = (z — H)" 1 and Gq(z) = (z — Ho) are the 
resolvents of H and Hq. T(z), V, G{z) and Go(z) are also related by 

G(z) = Go(z) + G (z)T(z)Go(z) (16) 
T(z)G (z) = VG(z) . (17) 

We shall see that the matrix elements of the resolvents in coordinate representation are singular on the real positive 
axis, where there is a branch cut, and at poles on the negative real axis (bound states). Further singularities may 
occur by analytical continuation on the second energy sheet. 

Note that the operators of scattering theory have abstract or parameterized versions 0] . Confusion may arise if they 
are not properly distinguished. The relation between abstract and parameterized operators is found by acting with 
|l3| ) on a square integrable state. The resulting infinite-time limits can be substituted by the following limits, see e.g. 

r° 

n+=\ime dte et e iHt/h e- lHat l n (18) 

poo 

n_=lime/ dte- £t e lHt/h e- lHot/h . (19) 
Integrating, and introducing a closure relation in momentum, 

/oo 
dpn(E p ±iO)\p)(p\ , (20) 
-oo 

/oo 
dpT(E p ±iO)\p){p\. (21) 
-oo 

The action of these operators on plane waves is now well defined. In particular, the improper eigenvectors of H are 
obtained by acting with the parameterized Moller operators on the plane waves, 

- Q(E P ± i0)\p) = \p) + 1 T(E p ± i0)\p) , (22) 

rjp ± lU — Hq 
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where E p = p 2 /(2m) is the energy of the plane wave and of the corresponding eigenstate of H. This is the Lippmann- 
Schwinger integral equation for the states \p ), which are composed by a "free" plane wave and a "scattering" wave. 
To evaluate the coordinate representation and the asymptotic behaviour of the states at large distances, the matrix 
elements of the free-motion resolvent are required, 

1 • < T ^L e ±i\ P \\x- x '\/n_ (23) 



x 1 E p ±iO- H Q < 1 ^h\p\ 

( p3| ) is obtained by introducing a resolution of unity in momentum representation and using contour integration in the 
complex momentum plane. Note that the two ways of approaching the real axis in (^3|), from below or from above, 
imply different boundary conditions at large \x\ for the two states in (p2|): the scattering wave of \p + ) is formed by 
outgoing plane waves moving off the potential region, whereas the scattering wave of \p~) involves incoming plane 
waves towards the potential region. 

Since the plane waves \p) form a complete set, the following resolutions of the operators T and S can be 
introduced, 



n± = / ^ |p ± > ( P | , (24) 

J — oo 

/oo 
dplp^ip], (25) 
-oo 

S= dp dp'\p')(p>-\p+)(p\ . (26) 



Strictly speaking, the operators in (|24j-g6|) are not identical to the ones in ( |10|jl2| , |l3| ) since the domain of the former 
includes plane waves. However, when acting on Hilbert space states they are equivalent so that, to avoid a clumsy 
notation, the same symbols will be used. A momentum representation is therefore allowed for these operators, which 
in general involves distributions (generalized functions such as Dirac's delta or Cauchy's principal part). 

For real potential functions V{x) the norm is conserved throughout the collision, (</>i n |</>i n ) = = (</w|0out) ■ 

This means that the Moller operators are isometric, i.e., 

nln ± = i op . (27) 

As a consequence, 

(p ± \p' ± )=5(p-p'). (28) 

In general the Moller operators are not unitary because the bound states are not in their range. Contrast this to the 
the operator S: it conserves the norm too, but it is unitary because it maps the whole Hilbert space onto the whole 
Hilbert space, 

55 f = S*S = l op . (29) 

The scattering states i\) with incoming and outgoing asymptotes move far away from the potential so they are or- 
thogonal to the bound states {| at large (positive or negative) times. Since the overlap amplitude = is 
independent of time, the space of bound states B is orthogonal to the scattering states, namely to the range of the 
Moller operators. We shall always assume that the ranges of the two Moller operators are equal to the subspace of 
scattering states 1Z, and that the whole Hilbert space is the direct sum of the subspaces spanned by scattering and 
bound states, TL = 1Z B. This assumption is known as asymptotic completeness, 



fl ± fll = l op - A = / dp \p±) ( P ± | . (30) 

J — DC 

In this expression the "unitary deficiency" A is the projector onto the subspace of bound states, 

A = ^|$,)(^|. (31) 

3 

Taking matrix elements in (p6|), the momentum representation of S is given by 

{p\S\p') = S(p - p') - 2m5{E p - E p ,)(p\T(E p + i0)\p') . (32) 
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The collision conserves the energy, which is, asymptotically, kinetic energy. That is why S commutes with Hq and its 
matrix elements are proportional to an energy delta function. It is quite useful to factor out this delta function to 
define an on-the-energy-shell S(E) matrix. Using 

6(p-p') = ^5(E p -E> p )6 pp , , (33) 

where 5 pp i is the Kronecker delta, 

and defining the matrix elements of S, S a p, by 

(p\S\p') = \p\m~ l S{E p - E p )S sign{p)sign{p , } (E p ) , (35) 

one finds 

S , 8 ig n (p) S i gI1 (p')(-Bp) = S pp —(p\T{E p + iO)\p') , \p\ = \p'\ . (36) 

The subscripts a, = ± in the matrix elements S a p, denote the two possible "channels", which correspond to positive 
(+) or negative (— ) momentum. A difference between the one dimensional scattering on the full line (— oo < x < oo) 
and radial scattering on the half line, (0 < r < oo) is that in the former, the S matrix is a unitary 2x2 matrix while 
in the later is a complex number of unit modulus. 

C. Symmetries 

1. Time reversal invariance. 

This symmetry holds for real potentials. It implies 

S a p — S-f]- a . (37) 



2. Parity. 

Frequently the potentials are symmetrical with respect to its central position. In that case, 

S a fj — S- a -p ■ (38) 
D. Eigenstates of H 

The eigenstates of H given by the Lippmann-Schwinger integral equations (|22] ) behave asymptotically as a combi- 
nation of two plane waves with positive and negative momenta. The factors multiplying these plane waves are the 
reflection and transmission amplitudes according to the following table for asymptotic, long-distance behavior (assume 
for the time being that p > 0) 

1 j exp(ipx/h) + R l (p) exp(~ipx/h), if x ~ — oo /„„% 
hUz \T l {p)exp(ipx/h), ifx-oo, [ ' 

T r (p)exp(-ipx/h), if x oo ^ 



h 1 / 2 \ exp(— ipx/h) + R r (p) exp(ipx/h), if x ~ oo. 
For potentials of finite range that vanish outside [a, b] these are in fact exact expressions for x < a and x > b 
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For p < 0, however, the states determined by 
wave (x\p), and (x\(—p)~), with outgoing plane wave 



If p > 0, the boundary conditions in ( |39| ) define the states (x\p + ) corresponding to an incoming plane wave from 
the left, (x\p), while the boundary conditions in (ffol ) define the states (x\(—p) + ) corresponding to an incoming plane 
wave from the right, (x\ —p). T{p) and R(p), with superscripts / or r for right or left incidence, are the transmission 
and reflection amplitudes. A wave packet peaked around a given \p + ) would be dominated by the plane wave \p) 
before the collision, whereas, after the collision, there would be two packets, one reflected and one transmitted with 
probabilities \R(p)\ 2 and |T(p)| 2 , dominated by | — p) and Mp) respectively, see e.g. ||. 

and (4Q) correspond, respectively, to (x\p~), with outgoing plane 
— p). A wave packet formed around would be close 
to a plane wave \p) only after the collision occurs. To form this peculiar outgoing state, the incoming asymptote 
must combine waves incident from both sides of the potential barrier. This may of course be difficult to implement 
in practice, but it does not preclude the usefulness of these states as basis functions, and in general for applications 
where some control or selection of the products of the collision is required. 

The previous discussion should make clear that T(p), for p < 0, is not a standard transmission amplitude, because 
it is not the amplitude of the transmitted plane wave of the state \p + ), p < 0. However, it analytically continues 
the standard transmission amplitude (T(p) for p > 0) onto the p < domain, so the term "transmission amplitude" 
will be used irrespective of the sign of p, even though the physical meaning is different for the two possible signs. 
According to our notational convention, positive arguments of the amplitudes correspond always to states \p + ), while 
negative momentum arguments correspond to \p~) states. 



E. Relation between scattering amplitudes and basic operators 

Comparing the asymptotic (large \x\) behaviour of the states in ( |39|) and ( |4p| ) with the asymptotic behaviour in 
(^), the amplitudes R(p) and T(p) can be related to on-the-energy-shell elements of the transition matrix. We shall 
work out one case in detail: the scattering part of (x\p + ) for p > and x — > oo is 

> 

dx'(x\G (E p + i0)\x')(x'\T(E p + i0)\p) 
2nmi e lpx / n 



h p 
2nmi 



-ip*'/K( x >\T(E p + i0)\p) 



P 



x\p)(p\T + \p) . (41) 



Adding the free wave, hr 1 l' 2 e %px l h ', and comparing with (p9|), there results T l (p) = 1 — 2iirm(p\T + \p) / p for p > 0. 
The rest of the cases can be worked out similarly (Because of time reversal invariance, {p\T±\p) = {— p\T± \ — p), and 
T r (p) — T l (p). Therefore the superscript for the transmission amplitude will be dropped hereafter): 

T(p) = 1 - ?^l(p\T sign{p) \p) , 
, . 2min , , _ , . 

R (P) = —(-P\%ign(p)\p) , 

R r (p) = -^y L ,<p\T S i^ P )\-P)- (42) 

Some useful relations follow from (f42|), 

[T(-p)]* = T(p), p real . (43) 

R r \-p)* =R r <\p), preal. (44) 



From (|36|) and (42), the S matrix is given by 



S(p) S S(25)=l$$ P>0. 



G 



It is quite useful to consider S as a (matrix) function of p. In simple applications we only use S(p) with p > 0j^| but in 
fact we may also define S(p) for p < or even for complex p in te rms of the analytical continuations of the amplitudes 
T(p), R r (p), and R l (p). This extension will be discussed in sec. II G . 



F. The diagonal Sd matrix 

The S matrix (^J) has been obtained from the momentum representation of S using plane waves incident form one 
side, | ±p), but other on-shell matrices may be defined in terms of a different basis formed by combinations of | ±p). 
Of particular interest is the set j = 0, 1, that provides a diagonal matrix 

■*w-(T«° w )- < 46 > 

Unitarity implies that \Sj\ = 1, so the matrix elements may be written in terms of real eigenphase shifts 8j, Sj = e 2lSj . 
The \uj) are not mixed by the collision; these incident states produce an outgoing combination equal to the incident 
one, except for a phase factor. The diagonal Sd matrix is most advantageous for parity invariant potentials, since the 
linear combinations become simply even and odd wavefunctions, 

\ U0 )=2-V 2 (\p) + \-p)), (47) 
\u 1 )=2-^(\p)-\-p}). (48) 

From the asymptotic behaviour of \Uj) = and | ±p + ) we may relate reflection and transmission amplitudes 

for even potentials to the eigenphase shifts, 

R{p) = 2- 1 {e 2lSa - e 2lSl ) , (49) 
Tip) = 2- 1 (e 2lSf > +e 2lSl ) . (50) 

(Eq. ( p0| ) is in fact valid for arbitrary potentials.) The boundary conditions for the states \u^) arc 

/2\ 1/2 

lim (x\u£) = e lSo - cos(-px/h + 5 ) , 



/2\ 1/2 

lim (x\ut) = e lSo - cos(px/h + S ) , 

V' 1 / 

/2\ 1/2 

lim (x\uf) = ie lSl - I sin(px/h - Si) , 
/2\ 1/2 

lim (x|«+) = ie lSl - sm(px/h + S x ) . (51) 

It will be convenient for later manipulations to drop the constant complex phase factors and define real cigcnfunctions 
of H as 

(a#o> =e-* s °(x\u Q +) , 

(ar|^) = -ie~' lSl (x\ Ul +) . (52) 



3 Keep in mind that p > in the arguments of S or of the scattering amplitudes does not mean 
According to the sign convention described in 
parts: 



[ID 



incidence from the right", 
it means that the amplitudes correspond to states with outgoing scattering 



for right incidence, and | —p) for left incidence. 
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G. Complex momentum 



The properties of T(p) as a function of the complex momentum p are of importance for many applications S. Let 
the potential function V(x) be such that 

p oo 

dx\V{x)\(l + x 2 ) < oo . (53) 

3 

Then T(p) is meromorphic in Imp > with a finite number of simple poles i/3%, i/?2, , iftni Pj > 0, on the 
imaginary axis. The numbers —0^ /(2m) are the eigenvalues of H. Moreover, 

T(p) = 1 + 0(l/p) as \p\ ->■ oo , Imp > , (54) 

and there can only be a zero at the real axis, at p = 0, 

\T(p)\ >0 lmp>0, p^O. (55) 

In the generic case T(0) = 0, and 

T(p) = 7p + o(p), 7 ^ 0, asp^ 0, Imp > . (56) 
Since T(p) is meromorphic and it does not have zeros in the upper plane, the integral 

1 f dlnTjp) 

7T- / d P ~ A = ~ n b ( 57 ) 

2tti J a dp 

along the contour A consisting of [— R, — e], [e, R], a semicircle of radius e around the origin, and a large semicircle 
of radius R in the upper half-plane, provides, according to a theorem of complex plane integration, the number of 
zeros (none in this case) minus the number of poles of T(j>) enclosed (the bound states). The integral may also be 
evaluated using fl43"|) , and jSq) ; this gives 2i$t(-R) — 2«$ T (e) — Mr, where $t(p) is the phase of T, 

T(p) = |T(p)|ex P (i$ T ). (58) 

Combining the two results, 

$ r (Q)-$ T (oo) = 7r(n 6 -l/2), (59) 

which is Levinson's theorem for the case T(p = 0) = 0. Otherwise, there is no — in contribution from the small 
semicircle and the phase difference becomes just nub- The convention followed is that <&t(oo) = 0, so the theorem 
establishes the value of <&t(0). 

The possibility to continue analytically T(p) to the lower half plane will depend of the potential considered |io| |. 
Here we shall assume that the continuation can be performed (this is the case for example for potentials of finite 
range) and discuss the properties that these continuations must obey. From T^(z) = T(z*) and the relations ( (42| ) we 
find 

(R r > l (p)T = R r ' l (-p*) , (60) 
(T(p))* = T(-p*) , (61) 

so that if there is a pole of T(p) in the fourth quadrant at pa — ipi {pr,pi > 0), there must be also a pole in the third 
quadrant at —pR—ipi- For an isolated pole, and if pi is small, the phase of T(p) along the positive real line will increase 
rapidly by tt. From ( |50"| ) we see that poles of T(p) are generally poles of So or of S\. Since |T(p)| = | cos(<5o — <Si)|, if 
the resonance eigenphase shift also jumps by tt, while the other one remains approximately constant, the transmission 
probability along the real axis will pass across a maximum (1), or a minimum (0), or both, depending on the initial 
phase difference of the two eigenphase shifts. The above simplified picture will be blurred if the resonances are very 
close to each other, or the pole is far from the real line. 
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H. Unitarity and its consequences 



The unitarity of the collision S matrix, SS^ = S'S = 1, reflects the conservation of norm in the collision. It provides 
two relations: from the diagonal elements, 

\T(p)f + \R- l (p)\ 2 = l, (62) 

and from non diagonal ones, 

T(p)[R l (p)}* + [T(p)]*R r (p)=Q, preal. (63) 
Eq. (|63| ) leads to a relation for the phases, 

2$ T - $ Rr - = (2n + 1)tt n = 0,±1,±2,... , (64) 

where, as in flsjj), 

iT«'(p) = \R r ' l (p)\e^n l W . (65) 



III. A MEASURE OF THE COLLISION DURATION: THE DWELL TIME 



In classical mechanics the quantity 

pt2 pb 

T D (a,b;t 1 ,t 2 )ciassicai = dt dxg(x,t), (66) 

J ti J a 

where g(x, t) is the probability density of an ensemble of independent particles, is the avera ge over the ensemble of the 
time that each particle trajectory spends between a and b within the time window [^1,^2] PI- ^ n other words, this is 
an average "dwell" or "sojourn" time in the selected space-time region.^ Its formal quantum mechanical counterpart 
is 

pti pb 

T D {a,b;t x ,h)^) = dt dx \ip(x, t)\ 2 . (67) 

Jt x J a 

In principle the coordinates a, b > a and the instants t\ and t% > t\ are arbitrary but most often a and b are chosen 
so that V(x) is zero or negligible for x < a and b > a. Hereafter t\ will be, by default, —00, or occasionally 0, an 
initial preparation time; and t 2 = 00. 

In spite of the formal similarity of the classical and quantum expressions, the interpretation of ( |67j ) as a "mean 
time" spent in the region [a, b], [ti, t%\ by quantum particles is not straightforward, since in the standard interpretation 
of the quantum mechanical formalism there are no trajectories and therefore there is no obvious way to assign a time 
(duration) of presence to a given member of the ensemble of particles associated with the quantum state. There 
are however several arguments that provide ( |67| ) by extending to the quantum case the classical dwell time, e.g. via 
Feynmann path integrals Jl3], causal or Bohm trajectories jl5]] , or as an expectation value of a hermitian sojourn time 
operator p6| , see also Chapter ... for an interpretation in terms of weak measurements. Irrespective of a hypothetical 
statistical interpretation of the dwell time in terms of individual members of the ensemble, the dwell time is at the 
very least a characteristic quantity of the ensemble represented by the state ifi, that quantifies the duration of the 
wave packet collision. For example, the dwell time is considered an important parameter in high speed applications 
of mesoscopic semiconductor structures Q. 

td can be written in several ways, in particular as 

/oo 
dtP ab {t) = (V>(t = 0)\T D \i/>(t = 0)) , (68) 
-00 



4 The concept of a "dwell time" for a finite space region in the stationary regime is due to Buttiker |12| . Previously, integrals 
of the form (|6q) had been used to define time delays by comparing the free motion to that with a scattering center and taking 
the limit of infinite volume, see e.g. 113] 
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where P a b{t) — J dxg(x,t). To is the sojourn time operator, 



T D = dte lHt/h D(a,b)e' lHt/n , (69) 



and D(a,b) is the projector onto the selected space region, 

D(a,b) = [ dx\x){x\ . (70) 

J a 

An experimental determination of the dwell time may be carried out by monitoring the time evolution of the probability 
inside the selected spatial region, p8[ . This is admitedly an indirect route, where the first moment of To, to, is 
obtained without having measured individual dwell times for the members of the ensemble. It remains to be seen if 
second and higher moments of To may be associated with some simple operational procedure. 

Let us now find other useful expressions for the dwell time. Integrating the continuity equation over x between a 
and b, and over time between — oo and t, P a b takes the form 

Pab(t)= f dt'[J(a,t')-J{b,t')}=[ dt! AJ{a,b,t') , (71) 

J — OO J — OO 

where J(x,£') is the current density, AJ(a,6) = J(a) — J(b), &nd the boundary condition P a b{— oo) = has been 
assumed. Substituting ( frjl ) into (|68|), one finds 

/OO pt poo pOO 

dt dt'AJ(t')= dt dt'H{t-t')f(t') (72) 
-oo J — oo J — oo J — oo 

= lim / dt'(t"-t')AJ(t') = lim t"P ab (t")-f dt't'AJ(t') . 

t"^aoJ_ 0o t"~,co J_ 00 



Unless P a b(t) decays faster than t , the dwell time wi ll div erge. The existence of a potential function leads generically 
to an asymptotic decay ~ t~ 3 , as discussed in section VC. However, for free motion the dwell time will diverge unless 
t he m omentum wave function vanishes at p — 0, because of the dependence ~ t 1 / 2 of the free motion propagator, see 
(141) below and the related discussion. In terms of the sojourn time operator (69) for H n , the possible divergence is 
due to a factor, 

To,h = £ r dpAp){p\D\ap){ap\ . (73) 

a =±J-oo \P\ 

In this and the following sections we shall limit ourselves in general to incoming asymptotes in the positive momentum 
channel (+) that vanish at p = 0, so that the dwell time for free motion does exist. This will allow to compare dwell 
times with and without potential and to define delay times. These states, with a bounded support in momentum 
space, have necessarily a Fourier transform in coordinate space that can only vanish at some set of points of measure 
zero. But this is not a problem since the total probability for positive positions tends to zero as t — > — oo, 

/•oo />0 

lim / dx\(x\4>- m (t))\ 2 = / dp|(p|0 in (O))| 2 (74) 

t — > — OO / / 

J a J — oo 

for any a and any <j)\ n 

Assuming that tP a b(t) — > as t — » oo, the dwell time ( |68| ) takes the local form 

/oo 
dt' [J(b,t') -J(a, t')]t' . (75) 
-oo 

Other expression for states incident in the positive momentum channel may be obtained by using resolutions of the 
identity in terms of the states |p + ), 

T D (a,b;^)= [ dp\{p\^ n (0))\ 2 ro(p) , (76) 
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where 



rdx\(x\ P +)\ 2 , 

r D {a,b;p) = ^ '/ I " , 77 

p/mti 

which suggests the interpretation of T£>(a, b;p) as a dwell time for particles of definite momentum p [|l2| . 

Suppose now that a < and b > are both far from the barrier region, before and after the barrier respectively, so 
that the first passage of the wave packet across a can be described accurately in terms of the free motion asymptote 
4>i n , while the passage of the transmitted and reflected wave packets can be evaluated with the asymptotic expressions 

Mb, t) = -±= jf dp (p|MO)) T (P) e^ b - Et ^ h , (78) 

Ma, t) = J™ dp (p|<MO)> R{p) er^+EDin (79) 

(For a potential with support between and d, b could be taken at the very barrier edge, b = d, but a cannot be 
0, because of the strong interference between the incident and reflected parts. \a\ should be much greater than the 
incident wave packet width in order to distinguish clearly the entrance passage from the reflected one.) Then, 

x: />oo 

dt'J T (b,t')= / dp\T(p)\ 2 \(p\fr n (0)}\ 2 = P T , 

oo JO 

dt'J^t 1 ) = / dp\{p\(p in {Q))\ 2 = 1 , 

oo J 



oo 



dt'J R (a,t') = - dp\R(p)\ 2 \(p\cl> in (0))\ 2 = -P R , (80) 



where the subscripts /, T and R in J/, Jt and Jr mean that </>;„, ipT and ipR have been used to calculate the fluxes. 
One can then write ( |7q ) as 

T D = p T (t)r-(ty: + p R (t)T\ (si) 

where 

t r dt'J T (b,t')t' 

/ t \out J—oo v ' / fQ<-)\ 

{)b = j T (b,e) ' (82) 

/oo 
dt'Ma^t', (83) 
-oo 

t _ - I-^dt' J R {a,t')t' 
{t)T = f?~dt>\J R {a,V)\ ■ iM) 

In each case the "average passage instant" is obtained by properly normalizing the fluxes. One may rightly wonder 
whether the notation and terminology used (as average passage times) are justified. The "averages" are taken over 
the current density J, a quantity that is not definite positive even for an incident wave packet without negative 
momentum components |20-22|. It turns out, however, that the above "averages" over J are equal to averages over a 



positively defined arrival time distribution (Kijowski's arrival time distribution) |23J, as will be discussed in Chapter 
... Models of detectors based on complex non-hermitian potentials lead also to these average times, delayed only by 
the small (dwell) time that the particle spends in the detector before being detected p3 . In the next section we shall 
relate these times to the "phase times" . 

Finally, note that (^) could be, and has been, used to partition the dwell time into transmission and reflection 
components [0,^5), see also the closely related aspproach of Olkhovsky and Recam i |2^| . The main drawback is 
that the defined entrance average instant is common for both contributions, see [|3,|8|[7which is not a correct in 



the classical ensemble limit, and may lead to negative transmission times 29 even in the classical case 
two-detector model avoided this problem by assigning different entrance instants for each member of the ensemble 
p8|| . The distinction between the dwell time and its components was first done by Buttiker fL2|| , and raised some 
controversy. As summarized in chapter 1, Muga, Brouard and Sala have emphasized the multiplicity of possible 
quantum partitionings versus the uniqueness of the classical case, and developed a systematic theory to generate 
partitionings with the correct classical limit. Some of these include interference terms that cannot be assigned to 
transmission or reflection but to both of them mAl . 
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IV. IMPORTANCE OF THE PHASES. TIME DELAYS. 



If the S matrix is known or simply one of the amplitudes R l or R r is given as a function of momentum and there 
are no bound states, necessary and sufficient conditions are known for a unique potential to exist, and there are well 
established construction procedures [^OP] . However a knowledge of the probabilities is not enough to determine the 
coefficients. The phases are associated with observable time dependent properties. 

Consider a wave-packet impinging from the left on a barrier potential located near x = (The exact barrier position 
is not important for our present purposes. Two typical choices for x = are the center of a symetrical barrier, or the 
left edge of a finite range potential). Let us take as before the spatial interval [a, 6] well outside the barrier, so that 
there is a clear separation between incoming and reflected passages. 

Since the incoming state is in the positive momentum channel, 

/>oo 

(a#i»(*)>=/ dp{x\p){p\^ n {Q))e^ x - Et ^\ (85) 
Jo 

and, applying the Moller operator fi +) 

/•oo 

<a#(i))=/ d P (x\p+)mn(0))e {ipx - Et)/h ■ (86) 



(If the zero of time is taken well before the wave packet interacts significantly with the barrier, the substitution 
(p\ip(0)) = (p|(/>in(0)) does not introduce any significant error.) 

Substituting (g5[), (|78|) and ( |79|) in the time averages (82-p4|), and using the standard expression for the current 
density, 

J(x,t) = -lJi>(x, y y^pll) , (87 ) 



dx 

the derivative of an energy Dirac's delta may be identified and then used to perform one of the momentum integrals. 
The results are 

(<>?"*= p- / dp\( P \^ n (0))\ 2 \T(p)\ 2 - lb-x +?i& T (p)] , (88) 

*T JO P 

<*>r = JT / dp\(p\^ n (0))\ 2 \R(p)\ 2 - [-a-x + K& R (p)} > (89) 

r R Jo P 



(t)a = / dp\(p\^ n (0))\ z -[a-x ) , (90) 
Jo P 

where the prime means derivative with respect to p, and 

ar = aM<&»(0)b>'<PlM0)>) • (91) 

These results do not require to assume a narrow packet in momentum representation. 
The quantity 



Ph 



(x ,b;p) = m[b-x +n^' T (p)]/p (92) 



in the integrand of ( p8[ ) consists of a time that a classical free particle with mass m and momentum p would spend 
from xq to b, plus the time delay mh& T (p)/p. Similarly, the term in brackets in 



-Ph 



(x ,a;p) = m [-a - x a + %& R (]p)] jp , (93) 



is a time spent by a classical particle that travels freely from xq to x — 0, where its momentum is instantly reversed, 
and from x = to a, plus a delay contribution. It is to be noted that unless a = —b the reference time associated with 
classical free motion is different in the transm ission and reflection cases. We shall see a consequence of this disparity 
afterwards when calculat ing average delays in IV B. 

Formally we may use ( |92|) and ( |93|) to define "phase times" for arbitrary values of a, b, and xq. In particular, for a 
finite range barrier between x = and d let us define 
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pt, , . md mh , , . , 
r£ h 0,d;p = — + —& T {p) , 94 
P P 

by substracting from T^ h (xo,d;p) the classical flight time between xq and 0, — mxo/p. These "extrapolated phase 
times" for traversal should not be overinterpreted as actual traversal times [^l|,^2| , not only because, as pointed out 
in Chapter 1, there is not a unique traversal time, but because a wave packet peaked around p is very broad in 
coordinate representation, so it is severely deformed before the hypothetical "entrance" instant t ent — \x Q \m/p, and 
at x = there is an important interference effect between incident and reflected components. The wavefunctions <j>i n 
and iJjr used to calculate the fluxes Jj and Jr do not faithfully represent the actual wave, so that the average instants 
(p9],p0|) loose their physical meaning as average detection times. 



A. The Hartman effect 



Relation ( |88| ) is suitable for examining the "Hartman effect" |3^,|34j,|31J,^G|,g5[ . Hartman |33|] studied the evolution 
of a wave-packet with momentum distribution centered around p c , colliding with a rectangular barrier of height 
Vo > p 2 /(2m), and width d. He found three regions according to the value of d. For large barrier widths (opaque 
barrier conditions), the stationary phase time associated with p, under the barrier, goes to a constant, T^ h (x , d;p) — 
2m/ (pn) — xqui/p, independent of d, where 

k= [2m(V -E)] 1/2 /h . (95) 

When transmission is dominated by momentum components below the barrier, the transmitted wave-packet seems 
to traverse the potential region in a time interval independent of d. This is the "Hartman effect" . If d is increased 
further, plane waves with momentum above the barrier height dominate the transmission, and classical behaviour 
results, i.e., time grows linearly with d. Finally, for small barrier widths, Hartman defined a "thin barrier region" 
where the phase time depends generally on d. 

To be more specific, let us consider the initial Gaussian wave-packet 



(z|<M0)) 



1 



2ttS 2 



1/4 



exp [ip c x/h - (a? - x c ) 2 /{A5 2 )} , (96) 



of average momentum p c = hk c and spatial width (square root of the variance) S. Here xq becomes equal the wave 
packet center x c . The initial momentum distribution is a Gaussian distribution with variance a 2 — [?i/(2<5)] 2 . We 
assume that p c >> a 2 so that the truncation at p = in ( |88| ) is not significant. For an energy distribution peaked 
around E c < Vq the following results can be drawn p5| : 

If K c d = ^/2m(Vo — E c )d/h >> 1, (t)^ ut does not vary appreciably when d increases, thus showing Hartman effect. 
When d is sufficiently large, the components of the wave-packet under the barrier are so strongly depressed by |T(p)| 2 
that higher momenta start to dominate, and (t)^ ut grows almost linearly, as one expects classically. As S is increased, 
larger values of d are needed to pass from the first regime to the second one. An estimation of the value of d which gives 
the transition between Hartman effect and quasiclassical behaviour can be obtained for each value of S by equating 
the factor |T(p)| 2 |(|0in(p))| 2 for p — p c and for p = p ri where p r is the momentum of the first resonance above the 
barrier. This leads to the relation 

_ hy/-]n\T(p c )\ hy/^d 
° ~ i i i 1> v*') 

\Pr-Pc\ \Pr-Pc\ 

between 8 and d, that clearly separates quantum and quasiclassical behaviour. Also, for fixed 5, the transition is 
sharper at larger 6 as a consequence of the narrower momentum distribution. 

We have already warned the reader against a naive overinterpretation of the extrapolated phase time T^ h (0,d;p), 
which becomes 2m/ (pn) for the barrier traversal in the Hartman effect, mainly because of the strong deformation of 
the broad incident wavepacket. We could try to avoid the interpretational pitfalls of this quantity and look instead 
at the time (t)^ ut for a wave packet which is initially localized near the edge of the barrier, and with a small spatial 
width compared to the barrier length d, to identify the entrance time and the preparation instant with a tolerable 
small uncertainty. However, Low and Mende speculated p5| , and then Delgado and Muga have shown p6| , that 
this localization leads to the dominance of over-the-barrier components. Similar conclusions are drawn from a two 
detector model, one before and one after the barrier, when the detector before the barrier localizes the particle into 
a small spatial width compared to d [E8J . 
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B. The lifetime and delay time matrices 



The four delay times corresponding to reflection and transmission for right and left incidence form the delay time 
matrix introduced by Eisenbud in his thesis |37|] , 



At a p = Re 



1 dS a 



S a p dE 



(98) 



The matrix element At a p is the delay time in the appearance of the peak outgoing signal in channel (3, after the 
injection of a pulse narrowly peaked in momentum in channel a. The "delay" may in fact become negative as discussed 
already. These delay times have been traditionally obtained by means of the "stationary phase argument" . Let us 
rewrite the transmitted wave function as 

/>oo 

(a# r (t)) = h ~ 1/2 / <W xp/fi -^ t/R+i * T Win(0))inP)l • (99) 
Jo 

If the initial state is narrowly peaked around po, the integral will be appreciably different from zero only if the phase 
of the exponential function is stationary near p — p$. This implies a "spatial delay" with respect to the free- motion 
wave packet, 



Air = h- 



dp 



P=Po 



and a corresponding "time delay" 



At++(;p ) 



Km d$T 



P=Po 



Po dp 

The time delays are also related to the on-the-energy-shell lifetime matrix of Smith |38| , 

dS(Ey 



Q(E)=ihS{E)- 



dE 



S is unitary, so Q is Hermitian. Thus the diagonal matrix elements of Q are real and take the form 



}aa — J]] \S a l3\ 2 Ataf; 



(100) 



(101) 



(102) 



(103) 







Since the particle has a probability IS^I 2 to emerge in the channel f3, Q aa is the average delay experienced by the 
particle injected in channel a. 

We shall now relate the Q matrix with the "wave packet lifetime" , defined as the difference between dwell times 
with and without potential p8|,p9|, 



(Q) = T D,ip - Td,, 



(104) 



As before, the incidence is in the positive momentum channel. td,i)> is given by (81) whereas the dwell time for free 
motion is 



Tn, 



(t) out 



(t)Z 



dp|W in (0))| 2 -[6-a] , 



where, similarly to (|9C|), 



(t) out 



b,4>i, 



dplipl^mf-ib-xo] 



(105) 



(106) 



Since, by hypothesis, (t)" 1 ^ — (t)™^. , (Q) takes the form 



(Q> = dt d X (\( X m))\ 2 ~ \(^in(t))\ 2 ) 

JO J a 

= - (t>$J + Pn[(t)Z% - (t)° b %J 



(107) 
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Substituting all the integral expressions obtained for the passage times, and writting c = 



—a — b, 



(Q) = nJ OO dp^\{p\4, i n(0))\ 2 [^'t\t(p)\ 2 + + \R(p)\ 2 } ■ (108) 

Note the term proportional to c in the reflection part. It arises because of the mismatch between the free motion 
reference times used to define the reflection and transmission time delays when c 7^ 0. Choosing c — 0, (Q) represents 
the weighted momentum average of the mean delay for each momentum,^] 

y* 00 

(Q) = / dp\(p\0 in (O))\ 2 Q(E) ++ . (109) 



The eigenvalues of Q have been used as good indicators of resonances J4lJ, see [V C below, and may be interpreted for 
symmetrical potentials as the delays associated with symmetrical or antisymmetrical bilateral incidence fl40[ . However 
their operational interpretation in terms of individual measurements is puzzling. An asymptotic measurement of the 
arrival time at b in the transmission side could be done in principle for one of the identically prepared systems 
represented by the wave packet. Because of the coordinate spread of the wavepacket, however, there is a large 
uncertainty in the time that the same particle enters the region [a, b]. If a detector is placed at a before the collision 
occurs, the entrance time can be determined, but in general either the particle is destroyed or its behaviour afterwards 
is modified by the measurement. We are thus faced with an intrinsic difficulty to measure individual delays. This 
means that, at variance with other quantum mechanical averages which are interpreted as averages of the eigenvalues 



measured for the individual members of the ensemble, the operational meaning of (109) does not require to assign 
a lifetime to a given particle. It depends on the average times defined in (R2m3), wl 



• rich are measurable, at least in 
principle, by the time-of-flight technique (Other operational procedure making use of particle absorption along the 
chosen interval has been described by Golub et al. [Q). This peculiarity of the delay time was already noted by 
Goldrich and Wigner |^| . A consequence is that the ordinary quantum fluctuations around the average value are 
not operationally meaningful. Instead, the relevant fluctuations refer to variations of the average values themselves, 
corresponding to S matrix (or Hamiltonian) ensembles 



The trace of (102) in the on-shell space is related to the change in density of states Ap{E) = Tr[<5(E — H) — 5(E — Ho)] 
which is a fundamental quantity to characterize the continuous spectrum [ flo] according to the "spectral theorem" 
(The three dimensional elastic and multichannel versions of the spectral theorem have been extensively discussed and 
proven rigorously |46|.) 

Ap{E) = -7r- 1 ImTr[G(E + iO) - G (E + i0)] 



Jew 



, d<S> T (E) 

.--^J. (110) 



The second equality (spectral theorem) follows from a result of Dashen, Ma and Bernstein |M . To obtain the final 



expression, use has been made of (62) and (|64J) ]48|, see |49J for an alternative derivation consisting in evaluating Ap 
for a finite system and then going to infinity. Note that the maxima of the trace of Q may be used to identify resonance 
energies and widths |5(J . For further relations between the density of states and the dwell time see |5l|-j53f . Chapter 
... discusses the concept of local density of states and its relation to the Larmor clock and transport properties. 

C. Breit- Wigner resonances 

The simplest model of resonance behaviour is the Breit- Wigner model for an isolated resonance, 

By imposing unitarity to S, and assumming that A and the resonance parameters Eq and T are independent of E, it 
follows that A = A^ and 



5 Additional oscillatory terms, see e.g. |S2|]4C|l , appear when the no-interference condition between the reflected and incident 
wave packets is not imposed. 
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A 2 =TA 



(112) 



This means that the matrix A factorizes as A a p — 7^7^, and that it is proportional to a projector matrix P = A/r 
with eigenvalues 1 and 0. Thus, the equation (|112|) takes the form 



r = £u 2 - (H3) 



The corresponding Q matrix may now be written as 

Q = Pq m , (114) 

with eigenvalues q m and zero, where 

nr 

Qm ~ (E - E ) 2 + T 2 /4 (115) 

is the maximum value allowed for a diagonal element of Q. The Breit-Wigner model for S and Q can be generalized 
in various ways, in particular to account for multiple overlapping resonances EJ. 



D. Negative delays 



In partial wave analysis of three dimensional collisions with spherical potentials, the time delay has been used 
mainly as a way to characterize resonance scattering. One of the standard definitions of a resonance is a jump by 
7r in the eigenphases of the S matrix. In one dimensional collisions the time delay has been also used frequently to 
characterize (non-resonant) tunnelling, where it may become negative. In fact the different delay signs associated with 
the two types of effects, resonances and tunnelling, are not independent. In 3D it was soon understood by Wigner H] 
that the increases and decreases of the phase should balance each other. Since Levinson's theorem imposes a fixed 
phase difference from p — to 00, there must be intervals of negative delay to compensate for the phase increases 
associated with the resonances. A similar analysis applies in ID to the transmission amplitude. In Figure 111 the phase 
of the transmission amplitude for a square barrier is shown versus p for different values of the barrier width d. 
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FIG. 1. Phase of the transmission amplitude versus momentum for a square barrier of "height" Vo = 5 and for three different 
widths, d = 1 (solid line), 2 (short dashed line), and 3 (long dashed line), m = 1. (all quantities in atomic units) 
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As d increases, the scatering resonances "above the barrier" p > po — (2mVo) 1 ^ 2 become more dense and are defined 
better because of the approach of the resonance poles in the fourth quadrant to the real axis. The corresponding 
increases of the phase are compensated by a more and more negative delay in the tunneling region. 

Negative delays also arise if a pole of T(p) crosses the real axis upwards, when varying the interaction strength, to 
become a loosely bound state in the positive imaginary axis. Levinson's theorem, see (|59|), imposes then a sudden 
jump in the phase $y(0) that must be compensated by a strong negative slope. This effect is more important near 
threshold, i.e., when the pole is very close to the real axis |5f|. Similar effects have been described for non-bound 
state poles in complex potential scattering |5fj . 

Wigner also found a bound for the negative (partial wave) delay time of a potential of finite radious. Whereas 
positive delays can be arbitrarily large, negative delays are restricted by "causality conditions" |57|] . Some back-of- 
the-envelope causality arguments may however be misleading. For example, assume a barrier of length d, and let a 
coincide with the left edge and b with the right edge. If the total time t^ h (0, d) is to be positive, the delay cannot be 
more negative than the reference free time, 

md , „. 

At++ > , 116 

P 

see e.g. p8|| . In fact this bound may be violated, in particular at low energy in the proximity of a loosely bound state. 
This should not surprise the reader after our repeated warnings against an overinterpretation of the extrapolated 
time (0, d). The flaw in the argument is the assumption of positivity of . Nevertheless, rigorous bounds have 
been established by Wigner himself and various authors in 3D collisions, see p9|,p7| f or review. In ID collisions the 
following bound holds for even potentials with finite support between —b and b pq|30f| : 

TYl I ! 

At++ > — < -2b H \sm(2pb/h + 2S ) - sm(2pb/h + 25i)l 

P { 4p 

This may be proven by using the even and odd eigenfunctions (x\%j}j) introduced in (|52|), in particular the fact that 
J^dxip 2 > 0. We start by calculating the logarithmic derivative of (x\?po) at x = b from the known expression for 
the outer region, see (|5l|), 



_ d(x\ip }/dx 
L b = 



(a#o 

Taking the derivative of L\, with respect to p, 



x—b 



-|tan(p6/n + 5 ) . (118) 



d5o f tl dLfr 2/ 7 /* , r \ , 1 _■ in/ 7 /-t , r l , X 



-T = - --r 1 cos 2 (pb/H + S ) + - sin[2(pb/h + S ] + - \ . (119) 
dp { p dp Zp n ) 

The first term on the right hand side may also be written as 

— [{x\4> a ) E {x\ijj ) x - (x\^ ){x\^ )e, x ](x = b) , (120) 

TO 

where the subscripts E and x are shorthand notation for the derivatives with respect to E and x. Repeating the same 
operations for x = — b one finds that 

[(a#oM2#o)x - (x\^o)(x\^ )e,x](x = b) 
= -[{x\ipo) E (x\ipo} x - (x\ipo)(x\i/j ) E<x \(x = -b) . (121) 

We shall now prove that this is a positive quantity. Taking the derivative of the stationary Schrodinger equation with 
respect to energy one obtains for real eigenfunctions of H the identity [ |38| 

h 2 d 

(xW = «a#>(x|V>£,* - (a#Ma#>,) , (122) 



so that, using (121) 
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b ^2 

dx {x\ip ) 2 = — {{x\i> ) E {x\ipo)x - (x\ipo)(x\ipo)E,x) {x = b) . (123) 
-b rn 



Carrying out similar manipulations for the odd wavefunction (x\ipi), (117) is found as a consequence of the positivity 
of the probability to find the particle in the barrier region. 

According to this bound the negative delay may be arbitrarily large for small enough momenta and may diverge at 
p = 0, as it occurs when a bound state appears when making the potential more atractive pq ]. For the square barrier, 
which does not have bound states, the time advancement of the Hartman effect is less important and it is actually 
bound by ( |116| ). Thus, whereas the experiments looking for anomalously large traversal velocities ( "superluminal 
effects") have been frequently based on evanescent conditions in square barriers (tunneling), square wells with the 
proper depth may in fact lead to much larger advancement effects at threshold energies. 



V. TIME DEPENDENCE OF SURVIVAL PROBABILITY: EXPONENTIAL DECAY AND DEVIATIONS 

The quantum mechanical decay of unstable states can be described in many different ways ]6l| , |62| ]. In many 
theoretical works the emphasis has been on justifying the approximately valid exponential decay law. A possible 
treatment for the survival amplitude A(t, ip) = (?/'(O)|'0(i)) decomposes the state tjj by the usual resolution into proper 
and improper eigenstates of the Hamiltonian H, corresponding to bound and continuum states. Even though it 
contains all the information, this is not convenient in general either for calculational purposes or for rationalizing the 
decay behaviour in a simple manner, except in favourable circumstances where the integral is easily approximated and 
parameterized, e.g. for isolated resonances and particular initial states. An ideal description would handle arbitrarily 
complex initial states and potentials in simple terms, and allow for an understanding of both the dominant exponential 
decay and the deviations from it. 

Much progress in this direction has been achieved by representating A(t, ip) as a discrete sum over resonant terms 
p3| , p4[ . The discretization allows a clear identification and separation of the physically dominant contributions, 
different terms being important for different time regimes. Here we follow the treatment presented in [ |64|j65| . 

The survival amplitude A(t, ip) = (ip(0)\ifj(t)} requires the diagonal matrix elements of the unitary evolution operator 
e -iHt/h^ when this operator is expressed in terms of the resolvent, A(t,ip) takes the form 

A{t^) = ^\e~ mt ' h ^) 

= 1^—1 dqqi^—^-m = ± [ dqe-^l h M{q) , (124) 
Zirm J c z — H Zir J c 

where z — q 2 /2m is a complex energy and the contour C goes from ™oo to +oo passing above all the singularities of 
the resolvent due to the spectrum of H (discrete poles for bound states and the natural boundary of the real axis for 
the continuum), and 

M(q) = ^(^|-^T7 |V> • (125) 
m z — H 

The survival probability is to be calculated as S(t,4>) = \A(t,ifi)\ 2 . 



A. Predicted time behaviour 



The function M(q) is evaluated in the upper half g-plane and then analytically continued into the lower half plane. 
Provided that the continuation exists, M(q) has in general a set of core singularities, depending only on the potential, 
plus possibly other structural state-dependent singularities. It is then useful to deform the original integration contour 
to the diagonal T> of the second and fourth quadrants of the g-plane. This provides both physical insight by identifying 
the most relevant time dependence (exponential decay) of the survival, and a calculational advantage for the remainder, 
since for t > the exponential e~ lzt / h = e~ lq t /( 2mh ) is a real Gaussian on this diagonal. 

Let us assume that a pole expansion of the form 

%) = EA ( 126 ) 

is possible (higher order poles can be treated in a similar fashion). Here fc = 1, 2, 3 ■ • ■ indexes the poles. On deforming 
the q integration from contour C to T>, the residues of the poles q^ crossed in the fourth quadrant on carrying out 
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this deformation provide contributions to A(t) that decay exponentially with time, whereas the residues are purely 
oscillatory for poles in the upper half plane (bound states), 

E k (t) = ake-^/^V = a k e~< , (127) 

where 

u = q/f, f= (l-i)v^fiT*) ■ (128) 

becomes real along the diagonal T>. Independently of providing or not providing a residue, all poles contribute because 
of the integral along the diagonal. Each pole contribution is expressed in terms of the w- function, see f66| , as 

Dk(t) = --^-sign(Imufe) iy[sign(Imufe)itfe] . (129) 

The exponential term may be added to this contribution to give the compact result JG6[ . 

A{t) = 5^[£7 fc (t) + D fc (t)] - ^a k w(-u k ) . (130) 

k k 

(It is understood that E k (t) = for poles in the lower half plane that have not been crossed when deforming the 
contour.) The second expression is very useful for studying the short time behaviour, but the first one has the 
advantadge of separating explicitly the exponential decay, E k , from the "correction" D k , which is given in terms of 
the known entire function w parameterized by the pole position and time. Numerical values and asymptotic properties 
of this function for small or large times are easy to calculate. 

The above treatment is easily extended for an M(q) that includes an entire function in addition to the pole 
expansion. This would add to the w-functions the integral along T> of the entire function times a real Gaussian. 



B. Short time behaviour 



The short time behaviour of the quantum survival probability is easily analyzed in terms of the above formalism, 
which allows to classify several possible non exponential dependences. 

Many authors have described a short time t 2 dependence of the decay probability Pdccay = 1—5 provided the 
mean energy and second energy moment of these states exist, see in particular the work related to the "quantum 
Zeno paradox" |57J. Less attention has been paid to the short time behaviour if these conditions are not fulfilled. A 
formal treatment and examples by Moshinsky and coworkers suggest a t 1 / 2 dependence of the decay probability at 
short times ]6qj69fl . We shall clarify how these two seemingly different claims can be compatible, and describe other 
possible dependences. 



The Taylor series (A7) of the w functions in (130) gives a series in powers of t 1 / 2 



[2- l q k {l-i){tlmh)V 2 ) n 



A(t, "^T^ m + T) ' (131) 

k n=0 y 2 > 

This suggests a short time t 1 / 2 dependence of the decay probability, as claimed by Moshinsky and coworkers |]68| , |69| . 
On the other hand, the formal series based on expanding the evolution operator, 



2 



provides a t dependence, 



o dcca y = ^((m 2 m - (mw 2 ) + ■■■■ (133) 



However, the expectation values of H and/or higher powers of H may not exist. Several behaviours are possible 
depending on the existence of these moments. The question of the physical realizability of Hilbert space states with 
infinite first or second energy moments is subject to debate |70[ ]. We shall leave this debate aside, and determine the 
possible implications on the short time behavior. 
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Consider the first two derivatives of A at time t = first from ( |132| ) and then by assuming a general short time 
dependence of the form A ~ 1 + 6£ c , where 6 and c are finite constants, 







dt 


t=o 


dA 2 




dt 2 





= _l { ^ m)=bct c-l 



t=0 



(134) 
(135) 



If the mean energy of the initial state does not exist, a i 1 / 2 dependence of the decay probability is possible, see 
examples in |6f| and [ |69| . 

If the mean energy is finite so that dS/dt\ t= o = 2Re(dA/dt\t=o) = 0, then c > 1. This rules out a i 1 / 2 dependence 
of A sin ce a t 1 / 2 dependence implies an infinite time derivative of A at t = 0. The corresponding coefficient for t 1 / 2 
in (131) must vanish by compensation between the different pole contributions. 

The second derivative is only finite at time zero if c > 2. This means that if the first energy moment exists but not 
the second, a dependence t c where 1 < c < 2 is possible for A (and for the decay probability), in p articular t?l 2 . (The 
coefficient for t 1 ^ 2 must also vanish in this case.) Otherwise, one can expect that the series (132) will be effective at 
short times for states with finite moments (ip\H n \i/)) leading to a t 2 behaviour. Examples where t 3 ^ 2 and t 2 dominate 
the short time behaviour of Pdecay are provided in p5| . 



C. Large time behaviour 



At first sight the asymptotic expansion of the w-function for t ~ oo in the correction term to the exponential decay 
suggests a long time dependence of the survival probability as t -1 , but in fact the ge neral behaviour is t~ 3 because 
of the cancellation of all the t~ x contributions. Due to the exponential e ~ lzt / n i n ( |138[ ) the large t behaviour is 
dominated by the region around the origin. The origin is actually a saddle point for the steepest descent path for this 
exponential factor that crosses the origin along the diagonal T> of the second and fourth quadrants. By introducing 



u and / variables as in (128) the exponential becomes e~ u and u remains real along the steepest descent path. 

The resolvent matrix element (ip\(z — ff) -1 |"0) which is defined for Imq > (first energy sheet) has to be analytically 
continued into the lower half g-plane (or second sheet of the complex z plane) to allow for this type of analysis, which 
will be valid in particular for finite range potentials. Provided that the analytically continued function is analytical 
at the origin it has a Taylor series expansion 

(ip\(z - H)' 1 ^) =a + axq + a 2 q 2 + ... (136) 



with coefficients Oj depending on rp. But because of the (odd) q factor in ( |125[ ), the first term, ao, does not contribute 
to the integral ( |138| ). The asymptotic formula for the survival amplitude comes therefore from the second term and 
takes the form 



1 — i ( mh 



\e-^ h W - -^-a.f / duu 2 e- 2 = J—L ai ^ . (137) 



3/2 



This formal result depends on the validity of (|136j), and on the assumption that no additional contributions due to the 
deformation of the contour are to be considered asymptotically. In general the analytically continued matrix elements 
of the resolvent will have poles in the lower half g-plane that may be crossed when deforming the contour, but these 
can only yield contributions that decay exponentially with time, so they are negligible at long times. 
A similar analysis may be performed for the propagator (no bound states) [ [71| 

{x\e- im ' h \x') = — [ dql(q)e- izt / h , (138) 
2tt J c 

I{q) = ^{x\^—\x'), (139) 
m z — rl 

substituting M(q) by I(q). Quite generally, I(q) vanishes at q = 0, and a i 3 / 2 dependence results. An exception is 
free motion on the full line, where 

M 1 \~>\ - ~ lm r l\x-x'\q/Ti 



z — Hq qh 



e l \ x - x , (140) 
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so that 1(0) — —i/h 7^ 0. As a consequence, the asymptotic behavi our o f the probability density for free motion on 
the full line is generically t . This is an important case in which ( 136 ) is not satisfied. Explicitly, by carrying out 
the integral in (138), the well known propagator 



(x\e~ lHot/h \x') = 



m \i/2 



(—) 

\ihtJ 



im(x—x') 2 /2ht 



(141) 



is obtained. A T 1 behaviour will also occur exceptionally when the potential allows for a zero energy pole of the 
resolvent. 

The free-motion probability density may decay faster than t _1 when the momentum amplitude (p\ip) vanishes at 
p = 0, so that the q~ x singularity is cancelled, see Figure 2. The exceptional cases of decay slower than i _1 has been 
studied by Unnikrishnan |J72| . 
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FIG. 2. dln\{x\ip(t))\ 2 /dlat versus dint for two different wave packets: one of them vanishes at p = 0, 
(p|^(0)) = C(l - e - ap2/n2 )e- s2{p - po)2/n2 - ipxo/n e{p) (solid line), and the other one is a a Gaussian wave packet, 
(p\tp(0)) = C'e~ s ( p ~ p °) / R -wxo/K (dashed line). C and C' axe normalization constants; the parameters are po = 1, Xo = —10, 
a — 0.5, 5 = 1, x = 0, m = 1 (all quantities in atomic units). Note the asymptotic dependences of the probability densities: 
t~ 3 and t -1 respectively. 



VI. OTHER CHARACTERISTIC TIMES OF WAVE PROPAGATION 



In the previous section we have seen how contour deformation techniques in the complex plane allow to single out 
contributions to the survival amplitude from resonance poles. In general the integral that provides the time dependent 
wave function may involve other critical points, "structural" poles, saddle points, or branch points, that determine the 
transient and asymptotic behaviour of the wave propagation. It is frequently possible to write explicit expressions or 
asymptotic expansions for the contributions of these critical points. In simple cases the effect of (the dominant term 
of) one of the critical points provides already a good approximation and a simple picture emerges, where characteristic 
times or velocities for the arrival of the main signal may be identified. Also typical is the transition from the dominance 
of one critical point to another, which may lead to a change in qualitative behaviour and to a characteristic time 
for the transition. The pioneering work in this direction is due to Stevens 73 741, who followed the techniques that 
Sommerfeld and Brillouin introduced in their study of the propagation of light in dispersive media |75fl . Examples of 
application to quantum scattering off a square barrier and a separable potential may be found in [[76|| and |f77f . Here 
we shall examine, following jTq] , the somewhat simplified case corresponding to a point source producing evanescent 
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waves following |7q| . This is not a "scattering problem" in the standard sense, but it illustrates quite clearly the 
techniques and concepts involved in more conventional scattering problems, and in other time dependent phenomena 
where a stationary state is achieved after a transient behaviour. 

In order to summarize essential aspects of the time dependence of wave phenomena other characteristic velocities 
or times have been traditionally defined. (We'll see that some of them coincide with times associated with critical 
points.) The phase velocity, uj/k, is the velocity of constant phase points in the stationary wave (assume k > for 
the time being) 

gifcr-iwt _ ( 142 ) 

The boundary conditions, the superposition principle and the dispersion relation lu = uu(k) between the frequency 
lu and the wavenumber k determine the time evolution of the waves in a given medium. If a "group" is formed by 
superposition of stationary waves around a particular to, it propagates with the group velocity diu/dk. In dispersive 
media (where lu depends on fc), the group velocity can be smaller (normal dispersion) or greater (anomalous dispersion) 
than the phase velocity. It was soon understood that these velocities could be both greater than c for the propagation 
of light; Sommerfeld and Brillouin |f75| , studying the fields that result from an input step function modulated signal in 
a single Lorentz resonance medium, introduced other useful velocities, such as the velocity of the very first wavefront 
(equal to c), or the signal velocity for the propagation of the main front of the wave. 

The above description is however problematic for evanescent waves, characterized by imaginary wavenumbers 
instead of the real wavenumbers of propagating waves. The role played by the imaginary part of the group velocity 
dio/dk and the possible definition of a signal velocity in the evanescent case have been much discussed. Assume that a 
source is placed at x = and emits with frequency luq from t — on. If luq is above the cutoff frequency of the medium 
(the one that makes k = 0) a somewhat distorted but recognizable front propagates with the velocity corresponding 
to uiq. For the dimensionless Schrodinger equation 

+ (143) 

the dispersion relation takes the form 

lu = 1 + k 2 , (144) 

and the signal propagation velocity for the main front is equal to the group velocity, v p = {dw / dk)^ = "2(ujq — l) 1 / 2 . 
In other words, at some distance x form the source, the amplitude behaves, in first approximation, as 

ip(x, t) sa e~ lulat e +lk ° x <d{t - xv p ) (145) 

where fco = (luo — l) 1 ^ 2 is the wavenumber related to luq by the dispersion relation, and Q is the Heaviside (step) 
function. In the evanescent case, u>q < 1, a preliminary analysis by Stevens J73|]74| , |79| ] , following the contour de- 
formation techniques used by Brillouin and Sommerfeld suggested that a main front, moving now with velocity 
v m = 2(1 — loq) 1 / 2 = Im(duj/dk)u , and attenuated exponentially by exp^cc) (where kq = (1 — luq) 1 ^ 2 ), could be also 
identified, 

i/>(x, t) a e - luJot e- KoX Q(t - xv m ) . (146) 
The contour for the integral defi ning the field evolution was deformed along the steepest descent path from the saddle 



point; and the main front (146) was associated with a residue due to the crossing of a pole at in^ by the steepest 
descent path. 

The result seemed to be supported by a different approximate analysis of Moretti based on the exact solution Jnj] , 
and by the fact that the time of arrival of the evanescent front, r = x/v m , had been found independently by Biittiker 
and Landauer J8^,[^] as a characteristic traversal time for tunnelling using rather different criteria (semiclassical 
arguments, the rotation of the electron spin in a weak magnetic field, and the transition from adiabatic to sudden 
regimes in an oscillating potential barrier). 

However, more accurate studies of the point source problem and other boundary conditions have shown that the 
contribution from the saddle point (due to frequency components above or at the frequency cutoff created by the 
sharp onset of the source emission), and possibly from other critical points (e.g. resonance poles when a square barrier 
is located in front of the source IfTql) are generally dominant at r, so tha t no sign of the luq— front is seen in the total 



wave density at that instant, see jSl|-)84 76 78y , and the subsection VIA below. 

Buttiker and Thomas reconsidered the signal sent out by a source which has a sharp onset in time p5fl . They 
proposed two approaches to enhance the monochromatic fronts compared to the forerunners due to the saddle. First, 
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the dominance of the high frequencies forerunners could be avoided if the source is frequency limited such that all 
frequencies of the source are within the evanescent case. Of course this makes the onset of the signal unsharp. A 
second option is not to limit the source but to frequency limit the detection. We can chose a detector that is tuned 
to the frequency of the source and that responds when the monochromatic front arrives. 

These two proposals and the sharp onset case were later implemented and examined in detail by Muga and Buttiker 
ifrsl . For a source with a sharp onset, they found that the traversal time t plays a basic and unexpected role in the 
transient regime. For strongly attenuating conditions (in the WKB-limit) the traversal time governs the appearance 
of the first main peak of the forerunner. In contrast, the transition from the forerunner to an asymptotic regime 
w hich is dominated by the monochromatic signal of the source is given by an exponentially long time, see more details 
in VIA below. If the source is frequency band limited such that it switches on gradually but still fast compared to 
the traversal time, the situation remains much the same as for the sharp source, except that now the transition fom 
the transient regime to the stationary regime occurs much faster, but still on an exponentially long time-scale. The 
situation changes if we permit the source to be switched on a time scale comparable to or larger than the traversal 
time for tunneling. Clearly, in this case a precise definition of the traversal time is not possible. But for such a source 
the transition from the transient regime to the asymptotic regime is now determined by the traversal time. Much the 
same picture emerges if we limit instead of the source the detector. Muga and Buttiker modell the detector response 
by means of a "spectrogram", a time- frequency representation of the wave function at a fixed point. As long as the 
frequency window of the detector is made sharp enough to determine the traversal time with accuracy, the detector 
response is dominated by the uppermost frequencies. In contrast, if the frequency window of the detector is made 
so narrow that the possible uncertainty in the determination of the traversal time is of the order of the traversal 
time itself, the detector sees a crossover from the transient regime to the monochromatic asymtotic regime at a time 
determined by the traversal time. 

Possibly, the fact that we can not determine the traversal time with an accuracy better than the traversal time 
itself tells us something fundamental about the tunneling time problem and is not a property of the two particular 
methods investigated. 



A. Role of the traversal time for a source with a sharp onset 



We shall obtain exact and approximate expr essions of the time dependent wave function for x > and t > 
corresponding to the Schrodinger equation ( 143 ) and the "source boundary condition" 



■4>(x = 0,t) =e~ luJot e(t) 



(147) 



in the evanescent case luq < 1. (A discussion of the physical meaning of "source boundary conditions" as compared 
to standard "initial value" conditions has been presented recently [p6f .) The solution may be constructed from its 
Fourier transform as 



dk 



1 



1 



Akx — ik 2 t 



(148) 



k + ino k — ino _ 

where the contour r + goes from ~oo to oo passing above the pole at ikq, and 

= (1-^ ) 1/2 ■ (149) 
The contour can be deformed along the steepest descent path from the saddle at k s = x/2t, the straight line 

ki = -k R +x/2t, (150) 

(kR and ki are the real and imaginary parts of k.) plus a small circle around the pole at in® after it has been crossed 
by the steepest descent path, for fixed x, at the critical time 



2k 

This procedure allows to recognize two w-functions J(3(| one for each integral, 

^(x,t) = le- it+ik ° t [w(-u' ) +«>(-<)] ■ 



(151) 



(152) 



Here. 
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2 l/2 * K 
1 + t 



(-<-?) 



(153) 



2V2 



1/2 



^0 ( « - T 



It is clear from the exact result (152,153), that r is an important parameter that appears naturally in the w-function 
arguments, and determines with kq the global properties of the solution. Its detailed role will be discussed next. 

The simplest approximation for ip{x,£) for times before r is to retain the dominant contribution of the saddle by 
putting k = k s in the denominators of (148 ) and integrating along the steepest descent path, 



^—it-\-ik t 



2i7T 1 /2 U 



1 



1 



(154) 



The average local instantaneous frequency for this saddle contribution is equal to the frequency of the saddle point 



LU s = l+ X 2 /U 2 . 

After the crossing of the pole in$ by the steepest descent path at t = r the residue 



—lUjQt —KQX 



9(t-r). 



has to be added to (|154j) 



ip(x, t) tp s (x, t) + tp (x, t) . 



(155) 



(156) 



(157) 



The solution given by Eq. (156) describes a monochromatic front which carries the signal into the evanescent medium. 
The conditions of validity of this approximation can be determined by examining the asymptotic series of the w(z) 
functions in (152) for large \z\, see the Appendix. In fact (157) is obtained from the dominant terms of these expansions. 
Large values of \z\ are obtained with large values of kq 7 t, or x, and also when t — > 0. Within the conditions that make 
the saddle approximation valid, the contribution of the pole is negligible. To see this more precisely let us examine 
the ratio between the modulus of the two contributions, 



R(t) 



2^/2 



-«o* f 3/2 (a ,2 /4f 2 + K 2y 



Its value at r is an exponentially small quantity, 

R(t = r) =e- KaX {2TTK x) 1/2 



(158) 



(159) 



In summary, fo r th e source with a sharp onset described here, the monochromatic front is not visible when the 
approximation (157) remains valid around t = r. A complementary analysis is carried out in Chapter .... 

However two very important observable features of the wave can be extracted easily from (157). The first one is 
the arrival of the transient front, characterized by its maximum density at tf = t/3 1 / 2 . This time is of the order of r, 
but the wave front that arrives does not oscillate with the pole frequency wo, but with the saddle point frequency u) s . 

The second observable feature that we can extract from (157) is the time scale for the attainment of the stationary 
regime, or equivalently, the duration t tr of the transient regime dominated by the saddle before the pole dominates. 
t tr can be identified formally as the time where the saddle and pole contributions are equal, R — 1. Because of (159) 
we shall assume r «t tr to obtain the explicit result 



ttr 



2 K 2^V2 



2/3 



(160) 



Finally, when xkq is small (~ 1), the saddle approximation describes correctly the very short time initial growth, but 
fails around r because the pole is within the width of the Gaussian centered at the saddle point. The pole cancels 
part of the Gaussian contribution so that the bump predicted by ip s at r/3 1 / 2 is not seen in this regime, r does not 
correspond to any sharply defined feature, but provides a valid rough estimate of the attainment of the stationary 
regime. 
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APPENDIX: APPENDIX: PROPERTIES OF ^-FUNCTIONS 



The w— function is an entire function defined in terms of the complementary error function as [ pq ] 

w(z) — eT z erfc(— iz) . 
w(z) is frequently recognized by its integral expression 



w(z) = — 



-du 



(Al) 



(A2) 



where T_ goes from — oo to oo passing below the pole at z. For Imz > this corresponds to an integral along the real 
axis. For Imz < the contribution of the residue has to be added, and for Imz = the integral becomes the principal 
part contribution along the real axis plus half the residue. From (A2) two important properties are deduced, 



and 



w(~z) — 2e z - w(z) 
w(z*) = [to(— z)]*. 



(A3) 
(A4) 



To obtain an asymptotic series as z — ► oo for Imz > one may expand (u — z) 1 around the origin (the radius of 
convergence is the distance from the origin to the pole, |z|) and integrate term by term. This provides 



w{z) 



1-3- 



(2m- 1) 



(2z 2 ) m 



Imz > 



(A5) 



which is a uniform expansion in the sector Imz > 0. For the sector Imz < (A3) gives 

I 3 • ... • (2m- 1)" 



w(z) 



(2z*y 



2e 



, Imz < . 



(A6) 



If z is in one of the bisectors then — z 2 is purely imaginary and the exponential becomes dominant. But right at 
the crossing of the real axis, Imz = 0, the exponential term is of order o(z~"), (all n), so that (A5) and (A6) are 
asymptotically equivalent as |z| — > oo. 
w(z) has the series expansion 



w 



{z)= f g 

w Vr(s + i)- 



(A7) 



The w-function is a particular case of the Moshinsky function jB8| , which can be regarded as "the basic propagator 
for a Schrodinger transient mode" [B7J. 
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